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ABSTRACT 

An angular correlation of low significance (2cr) is observed between 0-3 < z < 
0-5 QSOs and F < 23 galaxies. Overall, the cross-correlation function between 
82 intermediate-redshift (0-3 < z < 0-7), X-ray selected QSOs and y ^ 24 galax- 
ies is investigated, but no signal is detected for the z > 0-5 QSOs. After converting 
to an excess of galaxies physically associated with the QSO, this lack of strong cor- 
relation is shown to be consistent with the clustering of normal galaxies at the same 
moderate redshifts. Combined with previous observations, these results imply that the 
environments of radio-quiet QSOs do not undergo significant evolution with respect 
to the galaxy population over a wide range of redshifts (0 < z < 1-5). This is in 
marked contrast to the rapid increase in the richness of the environments associated 
with radio-loud QSOs over the same redshift range. 
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1 INTRODUCTION 

The study of QSO environments provides a direct measure 
of the bias in their distribution with respect to galaxies. This 
not only provides insights into the triggering/fuelling mech- 
anism of QSOs, but assists in the relating the results from 
QSO clustering studies (see e.g., Croom & Shanks 1996, La 
Franca, Andreani & Cristiani 1998) to large-scale structure 
studies carried out with galaxy redshift surveys. 

We have already investigated the galaxy cluster envi- 
ronments of low-redshift QSOs (Smith, Boyle & Maddox 
1995, hereinafter Paper I). The results of this study demon- 
strated that radio-quiet QSOs exist in average galaxy envi- 
ronments, consistent with these QSOs and galaxies having 
the same clustering amplitude. However, the observed ob- 
jects were all at low redshift {z < 0-3), and at low lumi- 
nosity {Mb > —24). In this paper we extend this work to 
higher redshifts and higher luminosities, probing the region 
of Mb,z space where the forthcoming large QSO surveys 
(2dF, Sloan) are targetted. 

Previous surveys have been limited by relatively small 
samples of radio-quiet objects, particularly in the redshift 
range 0-3 < z < 0-7. Yee and collaborators (Yee & Green 
1984, 1987, EUingson, Yee & Green 1991) have previously 
studied the environments of predominantly radio-loud QSOs 
at z < 0-6 and noted significant evolution in the environ- 
ments of QSOs towards rich clusters at higher redshift. For 
radio-quiet QSOs EUingson et al. (1991) measure an ampli- 
tude of IT ± 0-6 times the mean, zero-redshift galaxy au- 
tocorrelation amplitude and see no evidence for evolution 
over the redshift range, 0-25 < z < 0-6. At higher redshifts 
(1 < 2 < 1-5), most studies (Boyle & Couch 1993, Croom 



& Shanks 1998) find no significant clustering around radio- 
quiet QSOs. 

In comparison to earlier studies, the advantages of the 
present work are in the increased number of fields observed 
and the homogeneous. X-ray selection of QSOs. We are also 
assisted by having larger and fainter galaxy redshift surveys 
available with which to compare our results. 

In Section 2 we describe the data and its reduction. 
Section 3 summarises the techniques for measurement of the 
angular clustering amplitude and the estimation from this of 
a spatial clustering amplitude. Section 4 describes and dis- 
cusses the results and a summary is presented in Section 5. 

Except where explicitly stated to the contrary, an 
Einstein-de Sitter Universe has been assumed throughout. 
That is, flo — 1, qo =0-5 and Ao = 0. The Hubble constant 
is taken as 100/ikms~^ Mpc~^. 



2 DATA REDUCTION 
2.1 Observations 

We used deep l/-band CCD images to detect faint (V < 24) 
galaxies around QSOs in the redshift range 0-3 < z < 0-7. 
QSOs were drawn from the EMSS (Stocke et al. 1991) and 
CRSS (Boyle et al. 1997) serendipity surveys which use Ein- 
stein and ROSAT data respectively. Details of the EMSS 
data and the significance of X-ray selection are covered in 
Paper I. Here we simply note that X-ray selection ensures a 
fair sample of the general QSO population (Avni & Tanan- 
baum 1986). It is largely insensitive to bias in radio lumi- 
nosity and spectral index (Zamorani et al. 1981, Delia Ceca 
et al. 1994) and to moderate amounts of dust extinction 



2 R.J. Smith, B.J. Boyle and S.J. Maddox 



(Boyle & Di Matteo 1995) that can in principle seriously 
affect optical selection techniques. 

l^-band CCD images were obtained of 10-arcminute 
fields around 82 X-ray selected AGNs in the redshift range, 

0- 3 < z < 0-7 during two observing runs (17-23/09/95 & 15- 
21/04/96) at the 2-5- metre Isaac Newton Telescope (INT) 
on La Palma. The same 1024 x 1024 thinned Tek3 CCD cam- 
era was used on both occasions. The Tek3 chip has a peak 
quantum efficiency of 75 per cent at 6500A (>60 per cent 
throughout 4000-7700A) and mounted at the prime focus 
of the INT, a pixel scale of 0-5896 arcsec pixel" ^. The field- 
of-view corresponds to 2-l/i~^Mpc at z = 0-5. The CCD 
was read out in 'Quick' mode; read noise 5-9 electrons, gain 

1- 44 electrons ADU^^. Exposures varied from 15 minutes to 
Ihour depending on the QSO's redshift and seeing condi- 
tions at the time. In total 69 EMSS and 13 CRSS QSOs 
were observed. The redshift and absolute magnitude distri- 
butions of the sample are shown in Fig. 1. 



2.2 Source detection 

Raw images were bias subtracted, flat-fielded and coadded 
using the imred.ccdred package in IRAF. Source detection, 
extraction, classification and photometry were carried out 
using the SExtractor software (Bertin & Arnouts 1996). 
Each image was convolved with a Gaussian filter with full 
width at half maximum (fwhm) equal to the image's seeing, 
a background map determined and potential sources selected 
as connected groups of at least five pixels, each of which 
were individually greater than la above the background. 
For more details of the detection and extraction process, see 
Bertin & Arnouts (1996). 

Most fields contained saturated stars, artificial satellite 
trails and other defects. Such objects make detection of faint 
galaxies in their immediate vicinity difficult, so these regions 
were excised from the survey area in a semi-automated man- 
ner (see Smith 1998 for details). As an illustration, the holes 
used for field MS01084-I-383, a very low Galactic latitude 
{b — —24°) field significantly affected by bright stars, are 
shown in Fig. 2a. 

To optimise source detection and classification, it is nec- 
essary to know the seeing a priori. This information is used 
to choose a convolution filter for each frame and as the single 
input parameter into the star-galaxy separation algorithm 
(see below and Bertin & Arnouts 1996). The SExtractor 
software was therefore run on each frame twice. After the 
first run, the output catalogue contains the fwhm of every 
object in the frame allowing an accurate assessment of the 
seeing, which was taken as the modal fwhm value for un- 
saturated images. See Fig 2b. These seeing estimates were 
inspected individually. In a few cases with poor seeing the 
modal fwhm was not a robust estimate because of a broad 
scatter in the stellar locus. For these the seeing was esti- 
mated by inspection of the data frame and measuring in- 
dividual clearly stellar objects. In following this procedure 
it was consistently shown that more objects were detected 
in the second pass and the star-galaxy separation was more 
reliable. It was also confirmed that the second pass did not 
alter the assessment of the seeing. 
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Figure 1. (a) Redshift histogram of the 82 QSOs in the sam- 
ple, (b) Absolute magnitude histogram, (c) Redshift and abso- 
lute magnitude arc weakly correlated, a natural consequence of 
both evolution in the QSO luminosity function and this being a 
flux-limited sample. 



2.3 Photometric calibration 

Weather conditions during the observing runs did not merit 
calibration by photometric standard stars, so the CCD im- 
ages were calibrated by comparison with previously pub- 
lished, photometrically calibrated deep number-magnitude 
counts. We used the deep 1^-band galaxy counts of Small 
et al. (1995) and we took the stellar counts from the Galaxy 
model of Bahcall & Soneira (1980). 

By fitting the instrumental magnitude number counts 
from the CCD data to the star-plus-galaxy counts expected 
on the basis of the chosen models, the only free parame- 
ter was the photometric zero-point. The best fit was defined 
as that giving the minimum absolute deviation between the 
published and instrumental counts between V = 20 and 
0-5 magnitude brighter than the peak in the instrumental 
counts. The observed number counts were generally found 
to maintain an excellent fit to the model counts at V < 20 
until the CCD saturated at F w 18. The full catalogue was 
used in order to prevent errors in the star-galaxy separation 
from propagating into the calibration. One such calibration 
plot is shown in Fig. 2d. The observed number counts are 
shown as a histogram. The expected counts axe shown as 
the solid curve, with the stellar and galaxian counts plot- 
ted individually as dotted curves. Note that this is a low 
Galactic latitude (6 = —24°) field and the inclusion of the 
stellar counts becomes very important. For high Galactic 
latitude fields, the Small et al. galaxy counts alone provide 
an adequate fit to the shape of the observed counts. 

This method of calibration is subject to intrinsic varia- 
tions in the source number density, which are small at faint 
limits, and to the Poisson error on the counts. Assuming the 
bj < 24 galaxy auto-correlation of Roche et al. (1993) and 
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Figure 2. Four plots showing aspects of tlic reduction and photometric calibration process for one 
the top left: (a) All the sources detected in the field. The 1024 X 1024 pixel Tck CCD has been re 
purposes only. Large circles and squares show regions of the frame excised for any reason, usua 
Determination of the seeing parameter. The horizontal line shows the selected value for the stellar 
the star-galaxy classification flag against instrumental magnitude. A value of unity signifies a star 
these extremes lie on a non-linear probability function of being stellar. In all three plots, large 
final catalogue while small dots are objects which fall in regions of the frame excised for any 
significantly from a point-source profile (panel b) and arc successfully rejected (appear as small dotf^) 
(panel c). The final panel (d) shows determination of the photometric zero-point. The solid cur\ 
the histogram the observed counts. Calibration consists of sliding the histogram horizontally until 
dotted lines are the expected counts for stars and galaxies plotted separately. 



integrating over the CCD field of view, intrinsic galaxy clus- 
tering introduces an error of loss tiian 0-01 mag. Tlic star 
counts model is only claimed to have a 15 per cent accu- 
racy, but since stars typically constitute much less than five 
per cent of the sources, they also only introduce an error 
of order 0-01 mag. In the very low Galactic latitude fields, 
it may bo as much as 0-03 mag. With at least five hundred 
sources used in the calibration for each field, the Poisson 
error introduces an uncertainty of less than 0-05 mag. 



particular field, MS01084-f 383. From 
binned to 512 X 512 pixels for display 
proximity to a saturated star, (b) 
f^flun in pixels, (c ) A scatter plot of 
and^ro, a galaxl y. Values between 
points sniaW sources included in the 
)n. Note that=i|aturated stars deviate 
rather than miSBlassificd as galaxies 
e is the expected ^number counts and 
the best fitting position is found. The 



2.4 Star— galaxy separation 

The fields investigated covert wide range of Galactic lat 
itudes. Whilst the majoj ity have latitudes |&| > 40°, the: 
are a significant number nearer the Galactic plane and eight 
fields were observed at a h low as, \h\. < 30°., Stellar corita.mi 




nation of the galaxy calXguc ^jJLr,ft<t>(pa)^§,iffllBli}tHS 
any observed correlation (see Section 3.3). At low Galactic 
latitudes in particular, removal of the stars is therefore im- 
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portant. This may be done by excising the stars from the 
catalogue before estimating the correlation function, or by 
estimating the amount of contamination and re-scaling the 
correlation function to account for it. Both approaches were 
followed producing results entirely consistent within the sta- 
tistical errors. The results quoted herein are based on the 
procedure that used the star-gala^cy classification. 

The classification as stars or galaxies was achieved using 
the algorithm in the SExtractor package, which produces a 
classification flag with a value lying between (extended) 
and 1 (compact). A value of 0-8 was chosen as the boundary 
between stellar and galaxian images by inspection of plots 
such as Fig. 2c. This value includes the full width of the 
locus of points at 1 while excluding the cluster which fre- 
quently appears at ~ 0-5. Such a cluster is seen only weakly 
in the example shown here, but is much more pronounced 
for some fields. The algorithm's defaulting to values ~ 0-5 
is effectively a don't-know classification, but our cut of 0-8 
means such objects are called galaxies. 

Checks of the validity of the star-galaxy separation 
process, including the generation of artificial data using 
the artdata package in IRAF demonstrated that the star- 
galaxy classification procedure gave reliable results (see 
Smith 1998). 

2.5 The minimum observable angular scale 

Artificial data generated by the artdata package was also 
used to establish the effect of the bright QSO on detection 
of nearby faint galaxies. In order to do this 10 000 small 
(100 X 100 pixels) simulations were carried out using input 
parameters selected randomly from distributions defined by 
the complete CCD data set. The seeing, zero point, QSO 
magnitude and background level were all varied in this way 
(see Smith 1998 for details). Each image was furnished with 
sources down to one magnitude fainter than the CCD com- 
pleteness and then extracted and identified as usual. There 
was a marked failure to find galaxies within a few pixels of 
the modelled QSO. The ratio between the number of sim- 
ulated and detected sources was computed as a function of 
radius from the QSO. At distances greater than sixteen pix- 
els, the bright saturated object no longer had a detrimental 
effect and the source detection success rate was constant, 
but it decreased rapidly nearer the QSO. At less than five 
pixels from the QSO, the source detection rate was virtu- 
ally zero. No counts were therefore performed closer than 5 
pixels to the QSO. A best fitting, smooth polynomial was 
determined to give a correction factor for counts at separa- 
tions between 5 and 16 pixels. Though locally very signifi- 
cant, the 5 to 16 pixel anmilus is very small, so its overall 
impact on observed counts is also small and makes only a 
quantitative, not qualitative difference to the results. For ex- 
ample, counts out to one arcminute required multiplication 
by 1-0053 to compensate for the erroneously low counts in 
the central few arcseconds. 



3 ANALYSIS 

3.1 The correlation function 

The QSO-galaxy angular correlation function, Wqg{9), mea- 
sures the statistical excess of galaxies observed near the QSO 



over that expected for a random distribution of the same 
number of galaxies. In its standard form it may be defined 



n{e)5Vl = Ng[l^ wqg(e)]5^, 



(1) 



where n{d)6Q is the number of galaxies observed in solid 
angle 5f2 at angular separation 9 from the QSO. Ng is the 
background surface density of galaxies. The galaxy auto- 
correlation function may be defined similarly as the joint 

probability of finding a galaxy in both solid angle elements 

and SQ2 separated by angle 9. 



5P{9) =Ng[l+ Wgg{9)]SniSn2 



(2) 



The galaxy angular correlation function has been extensively 
studied (e.g., Groth & Peebles 1977, Maddox et al. 1990) and 
found to be well fit at scales less than 5° by a power-law of 
the form 



= Agg9 



1-7 



(3) 



where Agg is the angular covariance amplitude. The power- 
law slope is measured to be in the range 1-7 < 7 < 1-8. We 
here assume the QSO-galaxy correlation function to be of 
the same form. 



Wgg{9) 



■ Agg9'-' 



(4) 



We define Agg to be the amplitude at one degree except 
where explicitly stated to the contrary. 

Measurement of the angular covariance amplitude alone 
does not tell us about the physical environment of the QSO. 
We are integrating along the line of sight to the QSO, in- 
cluding information about the foreground and background 
galaxy populations. The more physically relevant spatial co- 
variance amplitude, Bqg , gives the strength of spatial clus- 
tering of galaxies actually associated with the QSO. The 
approach adopted here is to compare this with the galaxies' 
own clustering amplitude, Bgg, to determine the relative en- 
vironmental properties of QSOs and galaxies. 

By analogy to equation 1, the spatial cross-correlation 
function, ^qg{r), is defined as 



n{r)5V = pg[l + ^gg{r)]5V, 



(5) 



where n{r)5V are the counts in volume 5V at distance r from 
the QSO and pg is the average volume density of galaxies. 
The galaxy auto-correlation function is similarly defined as 



6P{9)=plll + ^gg{r)]dVi6V2 



(6) 



Limber's Equation (Limber 1953) relates the ampli- 
tudes of the angular and spatial auto-correlations and shows 
that equation 3 is a natural consequence of spatial clustering 
of the form 



Bggr~ 



(7) 



Measurements are quoted in the literature either as the spa- 
tial covariance amplitude Bgg, or as a correlation length 
ro.gg- QSO-galaxy cross-correlations have not been mea- 
sured to sufficient precision to allow a free fit, but are gen- 
erally assumed to follow the same simple power-law form, 
(e.g., Yee & Green 1987). 



£,gg{r) = BqgV ^ = 



(8) 
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Longair & Seldner (1979), hereinafter LS79, showed how 
to convert an observed Aqg to an implied Bqg . A similar ap- 
proach is adopted here with minor modifications. The con- 
version is performed in a statistical sense, since the galaxy 
redshifts arc not available. We therefore make several well- 
justified assumptions relating to the redshift distribution of 
gala^xies and the form and evolution of the galaxy clustering. 

First, all clusters are taken to be spherically symmetric. 
This may not bo so for individual cases, but is reasonable 
when averaged over a large sample. Secondly, we need to 
allow for the redshift distribution of galaxies, though not 
the actual redshift of any particular object. Previous ap- 
plications of this technique have predominantly derived an 
N{z) from a fixed luminosity function (LS79, Yates, Miller 
& Peacock 1989, Boyle & Couch 1993). Though it is possible 
to improve on this by using an evolving luminosity function 
(Ellis et al. 1996, Maddox et al. 1990), we have here adopted 
the parametric fit to observational redshift-survey data by 
Efstathiou (1995). (See also Baugh & Efstathiou 1993, Mad- 
dox, Efstathiou & Sutherland 1996.) Deep redshift surveys 
(e.g., Glazebrook et al. 1995) and studies of weak lensing 
around rich galaxy clusters (Small, Ellis & Fitchett 1994) 
mean the N{z) relation for B < 24 galaxies is now well de- 
termined, and this offers a more directly observational route 
to obtaining the required data. The parametric form for 
the galaxy number-redshift relation derived by Efstathiou 
(1995) is given by. 



diV 2 

— — oc z exp • 

dz 



Zc{bj) 



0-0113(6,7 - 17)^ ' 
0-0010(6,7 - 17)-^ 



3/2 " 



+ 0-0325 
0-0325 



17 < 6j < 22 

6,7 > 22 



(9a) 



(9b) 



The relation was normalised assuming a mean galaxy 6j — 
V = 0-5, from the same Small et al. (1995) V^-band counts 
as used for the photometric calibration in Section 2.3. 
The mean galaxy B — V colour was derived from the 
data of Fukugita, Shimasaku & Ichikawa (1995) and Driver 
et al. (1994), using the morphology mix shown in Table 1 
(see Prestage & Peacock 1988), and converted to bj — V via 
the colour equations of Blair & Gilmore (1982). The redshift 
distributions used are shown in Fig. 3. 

Finally LS79 assumed that the clustering is stable in 
proper coordinates. This is the case for virialized clusters 
where the structures behave as particles within the expand- 
ing Universe. On larger scales or for younger, non-virialized 
clusters, there is an intrinsic expansion or contraction within 
the cluster superimposed on the Hubble expansion. For a 
cluster simply expanding with the Hubble flow, the cluster- 
ing signal is constant in comoving coordinates. This is ex- 
pected in the case of sources that have a constant high bias, 
where the distribution of sources is determined by the initial 
conditions, largely independent of dynamical evolution. 

The rate of clustering evolution is here included as a 
free parameter, which makes the following changes to the 
derivation in LS79. Following Phillipps et al. (1978), we as- 
sume that only the amplitude varies with redshift, the form 
remaining constant. Equation 8 becomes, 

^qg{r,z) = Bo,qJ{z)r-''. (10) 

Let the typical cluster scale length, R, evolve as R{z) oc 



Table 1. Assumed galaxy morphology mix used in derivation of 
mean galaxy colours and K-corrections. 
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Figure 3. The galaxy redshift distributions adopted for the 
de-projection of the observed angular cross-correlation. They were 
derived from the functional form given by Efstathiou (1995). 
These may be compared with the QSO redshift distribution shown 
in Fig. 1. 



(1 -|- z)~^. From the definition of the correlation function 
(equation 1), the excess number of galaxies, iVe, in a shell 
around the QSO is 



N, = pg{z)6V{z)^qg{Riz),z). 



(11) 



This shell is then scaled with redshift according to the 
proposed evolution model. For observations to progres- 
sively higher redshifts, the volume sampled changes, making 
Pg{z) (X (1 + z)» and 5V oc R^ (x (1 + z)-^l^ . Thus, 



N,o,{l + zf{l + z)-''\f{z){l + z) 



(12) 



The first two terms arise from the change in volume sam- 
pled. The final term, which may be seen by substituting 
R{z) for r in equation 10, expresses the change in physi- 
cal scale sampled for the same comoving scale as a function 
of redshift. Assuming the number of galaxies physically in 
the cluster is a constant, iVe is a constant, and re-arranging, 
gives f{z) oc (1 -(- z)^''^~'^^~^ . Equation 10 is therefore re- 
written as 



(,qg{r,z) = Bo,qg{l + zYr~ 



: /3(3 - 7) - 3. 



(13) 



We substitute this form of the clustering evolution in 
place of that used by LS79, thus slightly modifying their 
equation expressing w{9) in terms of spatial parameters. The 
final expression relating angular and spatial amplitudes is 

then 



3-7/ 



{l + zy+~'^{mo,z). 



(14) 



The only difference between this and the LS79 equation is 
the power to which the redshift dependent term is raised. I-^ 
is a constant dependent only on the choice of slope for the 
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power law in equation 10 and D is the coordinate distance 
to the QSO. 



D : 



Ho{l + z) 



1 + 2; + VI + 2go-2 



1 + go2 + \/l + 2go^ 



(15) 



ip{mo, z) is an integral luminosity function, i.e., the number 
of galaxies per comoving cubic mega-parsec with apparent 
magnitude brighter than the flux limit of the sample, mo, 
at rcdshift, z. As described above, ip may be expressed as 
dN/AV , the ratio of an observed dN/Az and dV/dz, which 
is precisely specified by the chosen cosmology. 

The value of /3 then specifies the type of clustering evo- 
lution which is present. If/3=l,i^ = — 7 and the redshift 
dependence disappears giving comoving evolution. If /3 = 0, 
then 1/ = — 3 and the redshift term becomes (1 + 2;)"^"'"^, 
which is the case for stable clustering as in LS79. This evo- 
lution has also been parameterised in the literature in terms 
of the equivalent term e, where e = /9(7 — 3) and v = — (3+e) 
(e.g., Efstathiou et al. 1991). 

3.2 Practical estimators for Wqg{ff) and Aqg 

For a single field, the QSO-galaxy cross-correlation may be 

estimated by 



(16) 



where QG is the number of galaxies found in the range of 
separations (6* - M/2, 9 + Ae/2] from the QSO. QR is the 
rmmber of randomly scattered points in the same annulus. 
Random points were generated within the survey bound- 
aries taking full account of any holes or complicated survey 
geometry. In order to minimise random errors, many times 
as many random 'galaxies' were generated as there are real 
galaxies in the field, and the counts scaled appropriately. 

The clustering signal from any single QSO is very weak 
and has large associated uncertainties due to both Poisson 
errors and the intrinsic clustering of faint galaxies unrelated 
to the QSO. For most of the results presented here, the 
counts were coadded for all the fields within a given sub- 
sample (e.g., redshift slice), amounting to the pair- weighted 
average of j fields calculated individually, i.e.. 



«,,,(e) = - 1. 



(17) 



The angular covariancc amplitude was determined in 
two ways. If the full correlation function is plotted for a 
range of angular scales (e.g.. Fig. 4), the best fitting value 
was found for the assumed form (equation 4) of the corela- 
tion function. Alternatively, if counts are taken in a single 
broad annulus around the QSO, the assumed form may be 
integrated to give a scaling factor from the mean correlation 
in the annulus (the observed quantity) to the amplitude at 
any chosen scale. This is the approach adopted where Agg 
must be determined from a single QSO field, e.g.. Fig. 7. 
If QG and QR arc the galaxy and random counts in the 
range (Smim^max], then integrating equation 1 gives the ex- 
cess counts above the background 



QG-QR 



2TTdNgWggi9)de 



2neNgAgg ^ 6: 



(18a) 
(18b) 



using our standard assumed functional form for Wqg (9) . The 
background count, QR, is nNgiO^^y^ — &min)i which substi- 
tuted into the above equation gives. 



92 . 

■'max 
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-'mm 



q3-7 
max 



q3-7 
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(19) 



_3-7/QG_ 
2 \qR ^ 

If the angles are normalised such that 0max = 1, this formula 
gives the amplitude at the outer rim of the annulus. If they 
arc expressed in degrees, it would be the amplitude at one 

degree and so forth. 

3.3 Effects of poor star— galaxy separation 

The correlation function thus far defined is the true cor- 
relation function, wt{9). In practice the galaxy catalogue is 
likely to be contaminated by stars or incomplete through the 
mis-classification of galaxies. Since the background count is 
taken from the total number of galaxies in the field, both 
QG and QR are biased, which in turn biases the observed 
correlation function, Wo{9). 

In the case of stellar contamination leading to an over- 
estimate of the counts, it can straightforwardly be shown 
(e.g.. Smith 1998) that: 



(20) 



Mg and Ns are the total numbers of galaxies and stars re- 
spectively found in the field. Since this is a multiplicative 
factor, the power-law shape remains unchanged. Neglecting 
to perform any star galaxy separation, can therefore be cor- 
rected by rescaling the observed amplitude to its true value 
so long as the fractional contamination by stars is known. 

The alternative case is accidental exclusion of galaxies 
from the catalogue. Assuming any such failure is random, 
both QG and QR counts have the same fractional error and 
the observed correlation function remains unchanged. Fewer 
counts do of course imply a lower statistical significance on 
any measurement. We note that the random errors assump- 
tion may not be entirely precise. Failures of the star-galaxy 
separation are unlikely to be entirely uncorrclated, perhaps 
having greater effect on galaxies in particular environments, 
magnitude ranges or regions of the image. We have been 
unable to detect any such systematic errors in the current 
data. 

3.4 The Integral constraint, Ig 

The true background number density of galaxies is not 
known precisely, so was estimated from the data. The back- 
ground density was initially taken as the mean density 
within the particular field in question. If we then show that 
there is indeed a statistical excess of galaxies in the field, the 
estimated background value must have been biased high. 
The error may be retrospectively calculated from the ob- 
served correlation and the process iterated to self-consistent 
values. The fractional error in the background, Ib, or in- 
tegral constraint, causes the true background to be over- 
estimated by a factor 1 + Ib - Ib is determined by integrating 
the correlation function over the full area of the field used. 



Integrated excess counts / 27r^w(^)d^ 
Ib = = (21) 



True background 



Je„ 



2-Kede 
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assuming counts were taken from a circular region. For 
more complicated survey geometries, including holes in cov- 
erage, the quantity Ib is more easily obtained by numerical 
integration over each field individually. 

3.5 Errors on w{0) 

Estimation of the errors on a correlation function and quan- 
tities determined therefrom is a subject of active debate 
(e.g., Landy & Szalay 1993, Mo, Jing & Borner 1992 and 
references therein). The measured amplitude of the auto- 
correlation function at any given scale is dependent on the 
distribution of the galaxies over a wide range of separations. 
Any individual galaxy will exist as one member of many 
different pairs, so the pair counts in any one counting annu- 
lus are not independent of counts in the neighbouring annuli. 
As a result, despite essentially being a counting statistic, the 
error on the pair counts is not in general Poissonian. The im- 
portant distinction which must be drawn here is between an 
auto-correlation or full cross-correlation of two catalogues 
and the simplified geometry we are currently considering. 
This work simply aaimuthally averages the number density 
of galaxies as a function of the radial distance from the QSO 
position. As a result, the galaxies in all pairs arc independent 
so long as the fields do not overlap and the galaxy distribu- 
tion is fully defined by the QSO-galaxy cross-correlation. In 
reality galaxies are clustered, so an assumption that galax- 
ies' positions are independent is not true. It is however a 
very good approximation since any galaxy-galaxy cluster- 
ing is washed out around the counting annulus. The error 
on the galaxy counts should then be very close to the \/N 
Poisson error. This was confirmed by a full bootstrap anal- 
ysis (Barrow, Bhavsar & Sonoda 1984) of the errors (see 
Smith 1998) and so the Poisson error was adopted for this 
analysis. 

4 RESULTS 

4.1 Full data set, 0-3 < z < 0-7 

The cross-correlation, Wqg{0), for all the QSOs with V < 23 
galaxies is shown in Fig. 4. There is a positive correla- 
tion, though detected at a moderately weak level. At sub- 
arcminute separations there is a 2cr excess of galaxies over 

that expected for a random distribution and for the smallest 
counting annulus, approximately 22 arcsec, it is greater than 
3a. 

The V <23 cut was chosen to correspond to about one 
magnitude fainter than the characteristic absolute magni- 
tude. My — 5 log h, aX z ^ 0-4. It equals My — 5 log h at 
z ~ 0-55. This ensures that the peak of the galaxy N{z) re- 
lation lies in the same redshift range as the QSOs. My, as a 
function of redshift, was taken from the hj < 24 derived lu- 
minosity function of Ellis et al. (1996). Galaxy K-corrections 
(Oke & Sandage 1968) were taken from Pence (1976). Val- 
ues for diff'erent galaxy types were averaged according to the 
Table 1 morphology mix and a linear fit made over the range 
Q-3<z< 0-7. 

The predicted correlation functions for six difi^erent 
clustering amplitudes and evolution models have also been 
plotted in Fig. 4. Dashed lines are used for models stable in 
proper coordinates (/9 = 0, equation 13) and solid lines used 




-0.2 I 

12 3 4 

6 (arcmin) 

Figure 4. The cross-correlation function for all 82 0-3 < z < 0-7 
QSOs with y < 23 galaxies. Crosses are the observed val- 
ues. Open squares are modified by the integral constraint 
{Ig = 11-lAqg). The six curves show the expected signal for 
correlation lengths of ro = 10, 5-5 and 2/i~^Mpc respectively 
from top to bottom. The solid and dashed lines represent ro as a 
constant in comoving and proper coordinates. 

where ro is assumed to be stable in comoving coordinates 
(/3 = 1). To derive these curves, the correlation expected 
around each individual QSO was calculated (equation 14) 
assuming a power-law slope of 7 = 1-8 and the mean of the 
82 separate amplitudes (one for each QSO field) is plotted. 
In this way, projection efi^ects due to the redshift distribu- 
tion of the QSOs are accounted for. Whereas equation 14 is 
integrated over the generalised, mean galaxy N{z), this pro- 
cedure effectively integrates over the N{z) of the particular 
QSOs in the sample. The correlation lengths shown in Fig. 4 
are 10, 5-5 and 2/i~^Mpc, chosen to represent the following 
cases: 

ro — 2/!,~^Mpc: This is the correlation length of faint 
{B > 22) galaxies (e.g., Roche et al. 1993, Efstathiou 1995, 
Hudon & Lilly 1996). 

ro = 5-5/i~^Mpc: This is the present-day QSO-galaxy 
correlation length (Paper I). Also the galaxy correlation 
length, derived from bright, low-redshift galaxies has been 
measured by a variety of authors (e.g., Seldner & Peebles 
1978, Davis & Peebles 1983) to be in the range 3-4 to 
5-5h-^Mpc (Koo & Szalay 1984). 

ro — 10/i~^Mpc: This is an approximation to the clus- 
tering of galaxies around radio-loud QSOs. Yce & Green 
(1987) found correlation lengths of 8-7 and 17-5 /i"^Mpc 
for radio-loud QSOs in the ranges 0-3 < z < 0-5 and 
0-55 < z < 0-65 respectively. Using an enlarged data set, 
Ellingson et al. (1991) again found ro > 15/i^^Mpc for 
0-3 < 2 < 0-6, so 10/i~^Mpc may be viewed as a lower 
limit. 

For comparison. Loan, Wall & Lahav (1997) measured 
ro w 18 7i~^Mpc for a mixed sample of radio galaxies 
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and radio-loud QSOs. Results from the FIRST radio sur- 
vey axe giving the lower value of 8-12/i~^Mpc (Cress & 
Kamionkowsky 1998, Magliocchetti et al. 1998), possibly 
because the fainter flux limit is including many more star- 
burst objects than the Parkes and Greenbank surveys used 
in the Loan et al. paper. Dalton et al. (1994) measured 
ro w 14/i~^Mpc for the galaxy cluster-cluster correlation 
length, but if we imagine that all QSOs are found in galaxy 
clusters, a more appropriate comparison might be with the 
cluster-galaxy cross-correlation, measured by Lilje & Efs- 
tathiou (1988) to be ro = 8-8/i"^Mpc with 7 = 2-2. 

Prom previous studies, a general picture has emerged 
of Bgg being broadly similar to Bqg for radio-quiet QSOs. 
Ellingson et al. (1991) measured a QSO-galaxy correlation 
length of 5-6±3-2/i~^Mpc for their radio-quiet sample with 
no detectable variation over the range 0-3 < 2 < 0-6. Boyle 
& Couch (1993) detected no significant clustering around 
their optically selected QSOs (0-9 < z < 1-5). Though 
their result (ro = 0-9 ± 4-0 ^^"'Mpc) is consistent with a 
range of amplitudes up to and including the present-epoch 
galaxy covariance, it is inconsistent with the very rich en- 
vironments measured by Tyson (1986) for radio-loud QSOs 
in the same redshift range. Croom & Shanks (1998) actu- 
ally found a weak anti-correlation between bj < 23 galax- 
ies and optically selected QSOs. They show that this can 
be explained as an artifact caused by gravitational lensing, 
a conclusion consistent with the earlier study of Benitez & 
Martinez-Gonzalez (1997) which used an even brighter mag- 
nitude limit of bj < 20-5. 

In contrast, the results presented here are not consis- 
tent with the canonical ro = 5-5/i~^Mpc galaxy correla- 
tions. Also, they are greatly inconsistent with correlations 
as strong as have been claimed for radio-loud QSOs or 
the galaxy-cluster cross-correlation. We know, however that 
faint galaxies are not found in the same average environ- 
ments as bright galaxies either. Combining the results from 
five separate determinations of Ajg, Efstathiou (1995) found 
a constant comoving correlation length ro = 2/i~^Mpc to 
be a good fit for faint (22 < bj < 25-5) galaxies. A di- 
rect measurement of spatial clustering at moderate redshift 
from the Canada-France Redshift Survey (GFRS; Le Fevre 
et al. 1996) also gave the same value of ro = 2-0/i~^Mpc, 
assuming no evolution in comoving space between now and 
the mean redshift of the observations [z = 0-53). These clus- 
tering amplitudes are consistent with the observations pre- 
sented herein. As discussed by Lc Fevre et al. (1996), the 
inconsistency between observations at low and intermediate 
redshifts can be resolved either by assuming that a different 
population of galaxies is being probed by faint surveys, or 
conversely if the one population has undergone very signif- 
icant evolution over the range < 2: < 1. The Le Fevre 
et al. measurement of r{z = 0-53) — 1-33 ± 0-09 /i~^Mpc 
is consistent with ro = 5/i~^Mpc if the evolution param- 
eter, (3 = —1-6. This is a very rapid evolution. A value of 
/3 = —0-667 is predicted from f2o = 1 linear theory (Peebles 
1988). 

These results imply that radio-quiet QSOs apparently 
reside in environments consistent with those of the 'normal' 
galaxy population local to the QSO. We note however, that 
the observations are only 2a inconsistent with there being no 
correlation at all, i.e., ro = 0, and our detection of clustering 
is not as secure as that in large galaxy redshift surveys such 



Table 2. Covariance amplitudes and the corresponding correla- 
tion lengths for 0-3 < 2 < 0-7 QSOs and V < 23 galaxies under 
various assumed models of the clustering evolution, calculated 
from counts to one arcminute (equation 19). 



Covariance Amplitude 


Correlation length 


Inconsistency 






with 1 arcmin counts 


A (at 1 arcmin) 


ro (ft.~^Mpc) 




Observed 


Derived 




0-054 ± 0-020 


2-6 ±0-6 (f3 = 1) 




0-054 ± 0-020 


3-3 ± 0-6 {13 = 0) 




Derived 


Assumed 




0-606 (/3 = 1) 


10 


16 


0-398 (f3 = 0) 


10 


10 


0-207 {13 = 1) 


5-5 


5-0 


0-135 (/3 = 0) 


5-5 


2-7 


0-033 if3 = 1) 


2-0 


0-7 


0-022 (13 = 0) 


2-0 


1-1 


0-000 





1-8 



as the CFRS. Any 2D clustering results are equally consis- 
tent with an arbitrarily large correlation if we are prepared 
to accept the associated rapid evolution. 

The covariance amplitudes for 0-3 < z < 0-7 QSOs and 
V < 23 galaxies are summaxised in Table 2. (The actual 
observed counts may be found in Table 3.) In the first two 
rows, the observed amplitude is determined from counts at 
less than 1 arcmin and converted to the correlation length, 
ro- For the rest of table, assumed spatial-clustering mod- 
els are converted to Aqg via equation 14 for comparison. 
The final column gives the Poisson deviation required on the 
1 arcmin galaxy counts to produce the suggested covariance 
amplitude. The most probable correlation lengths for non or 
slowly evolving clustering are around 3 /i~^Mpc. Correlation 
lengths as great as ro — 5/i~^Mpc become acceptable only 
if there is very rapid evolution in the clustering amplitude. 
Clustering comparable to that reported for radio-loud QSOs 
is strongly ruled out. 



4.2 Subdivided data-sets 

In order to investigate the parameters in the evolution model 
and to reduce the effect of clustering dilution along the line 
of sight, we sub-divided the sample in redshift. We also re- 
peated the analysis for two sub-samples defined by absolute 
magnitude. 

Observed 1 arcmin pair counts for the low (0-3 < z < 
0-5) and high-redshift (0-5 < z < 0-7) samples are given in 
Table 3 and the cross-correlation functions are plotted in 
Fig. 5. Results from these samples demonstrate that most 
of the clustering seen in the full sample comes from the 
low-redshift fields, but that both sub-samples are consistent 
with either comoving or stable clustering with an ro of 2- 
3 ft"^Mpc. The 2 < 0-5 data are best fit by ro = 3-6 h"^Mpc, 
for virialized, stable clustering and ro — 2-9/i~^Mpc for co- 



moving clustering, assuming a 7 = 1-8 and the Aq 



0-0028 



least-squares determination of the amplitude. The slightly 
weaker Aqg = 0-0026 calculated from the one arcminute 
counts would give correlation lengths shorter than those 
quoted above by 0-1/i^^Mpc. The galaxy magnitude limits 
again correspond to M* — 5 log h+1 at the median redshift of 
the sub-sample. Inspection of Fig. 3 shows that in both cases 
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Figure 5. Cross correlations for z < 0-5 and z > 0-5 data plot- 
ted separately. In the upper, low-redshift panel, crosses are the 
raw, observed correlation and open squares are corrected for the 

integral constraint {11-0 Aqg). No correlation is detected in the 
highcr-redsliift sample, so no integral constraint is applicable. The 
dot-dashed curve in the upper panel is the measured amplitude 
oiAqg = 0026± 0-001, which gives ro = 2-8ft"iMpe (/3 = 1). A 
least-squares best fit to the observations gave the consistent value, 
Aqg = 0-0028. Compajrison curves are for 5-5 and 2 h~^Mpc as in 
Fig. 4. 

Table 3. Galaxy and random point raw counts at less than 
one arcminute from the QSO. cr is the Poisson significance of 
QG's difference from QR. Note that in two cases this is is a 
very small deficit rather than an excess, formally signifying an 
anti-correlation. After inclusion of the correction for 'detection 
success ratio' described in Section 2.5, these become a small, but 
insignificant positive correlation. 



QSO sub-sample 


Galaxy magnitudes 


QG 


QR 




0-3 < 2 < 0-7 


17 < y < 23 


1008 


948 


1-9 


0-3 < 2 < 0-5 


17 < y < 23 


691 


640 


1-9 


0-5 < 2 < 0-7 


17 < y < 24 


725 


729 


0-15 


My > -22 


17 < y < 23 


607 


548 


2-4 


My < -22 


17 < y < 23 


401 


403 


0-08 



the N{z) distributions peak within an appropriate redshift 
range. 

We can also combine the correlation length derived for 
the galaxy environments of z < 0-3 QSOs in Paper I with 
those derived here to give us a largo baseline over which 
to study the evolution of the clustering. A single measure- 
ment of the correlation defines a curve in the ro,f3 plane 



Figure 6. Values of ro,/3 which would be consistent with the 
present and previous observations. Shaded regions show la limits 
on our observations of Aqg . The solid line represents ^gg from the 
CFRS sample and the dashed line w{0)gg from the APM Galaxy 
Survey. 

giving all possible value pairs consistent with the obser- 
vation. Fig. 6 shows the regions of this plane bounded by 
the la uncertainties on our clustering amplitude determi- 
nations, < 2; < 0-3, 0-3 < z < 0-5 and 0-5 < 2 < 0-7. 
The high- redshift sample only defines an upper limit. Also 
plotted are two curves taken as examples of measure- 
ments. The dashed line is for bj < 20 galaxies from the APM 
Galaxy Survey (Baugh 1996) which have a median redshift, 
Zm = 0-13. The sohd line is for the CFRS (/ < 22-5), which 
has median redshift, Zm = 0-53. These are not necessar- 
ily absolutely comparable to our amplitudes, since although 
they cover very similar redshift ranges, they have difi^erent 
magnitude limits and used slightly flatter power-laws to flt 
the clustering amplitude. They are however in reasonably 
good agreement. 

The various observations may thus be reconciled if there 
has been moderately strong evolution in the clustering am- 
plitude, /3 < — 1 for both ^qg and ^gg. Such a conclusion must 
assume that the same single population is being sampled in 
each experiment. 

Alternatively, rather than summing the data over a pre- 
defined redshift range, in Fig. 7 the observed amplitude at 
one arcminute is plotted as a function of redshift for each 
QSO individually. The top panel shows the amplitude for 
each of the 82 fields with their individual error bars and 
three constant comoving correlation length models. In the 
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Figure 7. Observed QSO-galaxy cross-correlation amplitude at 
one arcminute using V < 23 galaxies. The two panels show the 
same data with an axis re-scaled for clarity. In the top panel, 
curves representing 10, 5-5 and 2/i~^Mpc comoving correlation 
functions are shown. In the lower panel, only 5-5 and 2/x~^Mpc 
are visible. Solid curves are for constant comoving evolution and 
dashed for constant proper correlation length. Filled boxes are 
each the average of eleven fields and the error bars give the one 
sigma standard deviation in each subsample. 

lower panel, the Aqg axis is re-scaled for clarity. Filled boxes 
are each the average of eleven fields, thereby adaptively sam- 
pling the redshift range according to the availability of data. 
The ro = 5-5/i~^Mpc and 2/i~^Mpc, comoving and stable 
clustering models are also plotted in the lower panel. This 
demonstrates that the summed data in the redshift slices are 
representative and are not strongly dependent on any single, 
or few, fields. There are significantly more positive points 
than negative in the z < 0-5, V < 23 data. When a similar 
plot is produced with a limiting galaxy magnitude V < 24, 
which we consider more appropriate for 0-5 < z < 0-7 QSOs, 
the data points are almost exactly evenly split around the 
Aqg — axis at all redshifts, giving agreement with the zero 
correlation seen in the lower panel of Fig. 5. 

The plot also shows that the predicted change in ampli- 
tude over the redshift range studied is small and the curves 
for different values of the parameter l3 are similar in shape. 
An order of magnitude increase in the number of QSO fields 
observed would be required to reliably distinguish between 
these two clustering evolution schemes. 

Finally, the clustering amplitude as a function of QSO 
absolute magnitude is shown in Fig. 8. Amplitudes are here 
plotted in terms of the spatial covariance to compensate for 



the different projection effects on QSOs of the same absolute 
magnitude, but different redshifts. The amplitude appears to 
decrease with increasing QSO luminosity. This is markedly 
different to the rapid increase in cluster richness with QSO 
luminosity seen by Yee & Green (1987) for radio- loud QSOs. 
Caution must be exercised when assigning any significance 
to the apparent variation seen in Fig. 8. Figs 7 and 8 are 
not independent since luminosity and redshift are related 
via the QSO luminosity evolution. Indeed, from Fig. 1 it is 
clear that 80 per cent of the My — 51og/i > —22 QSOs are 
also in the z < 0-5 sub-sample. What appears to be a corre- 
lation with luminosity may simply be the selection function 
being influenced by the evolving QSO luminosity function. 
Allowing for this, there is no evidence in the current data 
for a dependence on luminosity which is independent from 
that with redshift. 

4.3 Comparison with QSO and galaxy 
autocorrelations 

In the previous section we showed that, at intermediate red- 
shifts, 0-3 < z < 0-7, the inferred clustering scale length 
for faint galaxies around QSOs, ro = 2/i~^Mpc was sim- 
ilar to that derived from the auto-correlation function of 
faint galaxies, and thus, in the simplest interpretation, the 
QSOs and faint galaxies might exhibit the same bias. In con- 
trast, the most recent measurements of the QSO correlation 
function are consistent with a constant comoving amplitude 
ro = 5/i~^Mpc at z < 2-2 (Groom & Shanks 1996), suggest- 
ing QSOs are strongly biased with respect to the galaxy 
population if the galaxy scale length is ro = 2/i~^Mpc. 
However, the measurement of the QSO correlation length 
is dominated by QSOs with 1 < 2 < 2, with relatively few 
QSOs at 0-3 < 2 < 0-7 and certainly too few to rule out 
ro = 2/i~^Mpc for the QSO correlation length at these red- 
shifts. 

If the galaxy correlation length were to increase again 
at z > 1 then the QSO and galaxy distribution could be 
reconciled with the same value for the bias parameter (al- 
beit redshift dependent). There is some evidence that the 
galaxy correlation must increase again at higher redshift; 
angular correlations for HDF galaxies with photometric red- 
shifts shown a rapid increase at z ~ 2 (Magliocchetti & Mad- 
dox 1999, Arnouts et al 1999), and Giavalisco et al. (1998) 
measure ro — 6/i~^Mpc for galaxies at z ~ 3. Alternatively, 
the strong clustering seen for both QSOs and the Giavalisco 
et al. galaxies may simply be highlighting the danger of gen- 
eralising from any one type of galaxy to another when we 
do not understand galaxy formation and how the classes are 
inter-related. Such discrepancies are seen also for 'normal' 
galaxies, even within a single, optically selected survey. For 
example, both the luminosity functions (Lilly et al. 1995) 
and correlation functions (Le Fevre et al. 1996) of red and 
blue galaxies in the GFRS show significantly different evolu- 
tion over the redshift range < z < 0-8 and in comparison 
to the present epoch determinations. 

4.4 Future work 

It is clear that we do not yet understand the cluster environ- 
ments of QSOs with the same precision as we do for galaxies. 
New wide-field GCD cameras on large-aperture telescopes 
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Figure 8. QSO— galaxy spatial covariance amplitude, assuming 
/3 = 1, as a function of QSO absolute magnitude. Crosses are 
2 < 0-5, open boxes z > 0-5. Horizontal bars on the left of the plot 
indicate the amplitudes equivalent to correlation lengths, ro = 2, 
5-5 and 10/i~^Mpc. My > —22 is almost a subset of z < 0-5, so 
trends in Bqg from this plot and Fig. 7 are not independent. 



will be able to make some contribution to these studies, 
but the broad range in redshift of galaxies in a simply flux- 
limited sample places a significant limit on the sensitivity 
of this approach. The ever improving availability of redshift 
data f or faint galmcioa will go oomo way to allowing detailed 



studied of QSO pn-irirnnmpnts, though pvpn the iiprnming 
generation of sky surveys will not probe to sufficiently faint 
limits to obtain redshifts for typical M* galaxies at the 
redshifts investigated here. Photometric redshifts obtained 
from wide-field imaging in several pass-bands may prove the 
most efficient method to further the work described here. 
The potential for photometric redshifts at 2; < 2 has been 
proved through recent interest in the technique spurred by 
the Hubble Deep Field (e.g., Fernandez- Soto, Lanzetta & 
Yahfi 1998). 



5 SUMMARY 

Clustering of < 23 galaxies around X-ray-selected, 0-3 < 
z < 0-7 QSOs has been detected with a >2a signifi- 
cance. The amplitude of the correlation function is small, 
Aqg{l°) = 0-0020 ± 0-0008, implying a spatial clustering of 
amplitude ro « 3 h~^Mpc. Though not a secure detection of 
clustering around X-ray selected QSOs, it is a highly signif- 
icant rejection of rich cluster environments similar to those 
found for radio loud QSOs (Yee & Green 1987), and of the 
cluster environment around present day normal galaxies. 

When the sample was subdivided ai z — 0-5, the corre- 
lation was found to be stronger in the 2: < 0-5 data. No signif- 
icant excess of QSO-galaxy pairs over the mean background 
counts was observed at z > 0-5. The best-fitting amplitude 
for the low-redshift sample was ^,9(1°) = 0-0026 ± 0-001, 
which corresponds to ro « 3 ft^^Mpc if we assume there has 
been little evolution in the clustering since z = 0-5. 

These measured correlation amplitudes are in good 
agreement with measurements of the faint-galaxy correla- 
tion function (e.g., Roche et al. 1993, Le Fevre et al. 1996). 
Taken in conjunction with the results of Paper I, this implies 
that QSOs are likely to be populating average galaxy envi- 
ronments at all redshifts 2 < 0-5. Though no clustering was 
detected for the z > 0-5 subsample or in the high-redshift 
sample of Boyle & Couch (1993), these null results too are 
consistent with the same conclusion. This is in very marked 



contrast to the rapid evolution in the richness of galaxy clus- 
ters associated with radio-loud QSOs, which by z ~ 0-6 
are typically found in Abell richness class 1 clusters (Yee & 
Green 1984, EUingson et al. 1991). While further work re- 
mains to be done, the results here indicate that radio-quiet 
QSOs may not be strongly biased with respect to the galaxy 
population. This has importaint, and potentially positive, 
consequences for studies of large-scale structure from the 
forthcoming large QSO redshift surveys. 
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